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FACILITY: MATH LIBRARY 


: ABSTRACT 
This module contains routine MTH$CSQRT ~- compute COMPLEX square root. 


; VERSION: 0 
; HISTORY: 


; AUTHOR: : 
Jonathan M. Taylor, 20-JUL-77: Version 0 
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-SBTTL DECLARATIONS 


; : INCLUDE FILES: 


: EXTERNAL SYMBOLS: 


SA AAAAAAAOAO 
MEW" OOONOUS WO 


.GLOBL MTHSCABS 
0 *GLOBL MTHSSQRT_RS 
fi 3 
72 ; MACROS: 
0 73; NONE 
000 7 
800 «78 |; PSECT DECLARATIONS 
00000000 : -PSECT _MTHSCODE PIC, SHR, LONG, EXE, NOWRT 
0000 78 
4 79 ; 
000 80 ; EQUATED SYMBOLS: 
00000004 44 8 argadr = 4 ; offset from AP of arg address 
4 : 3 
00 84 ; OWN STORAGE: 
0000 85 ; NONE 
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MOOS MTHSCSQRT = compute COMPLEX square root ~$eE =} 382 94:59:38 ERTHRTL. SREIMTNCSORT MAR: 1 — 
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p++ 
; FUNCTIONAL DESCRIPTION: 


The square root of a complex number (r, i) is computed 
as follows: 


ROOT = SQRT(C(ABS(r) + CABS((r, 1))) / 2) 
Q= i / (2*ROOT) 


ue 
i 
000 9 
$8 3 ; r i CSQRT((r, i)) 
000 10 ; 
4 1 3 >= an (ROOT, Q) 
00 103; <0 = (Q, ROOT) 
909 <0 <0 (-@, -ROOT) 
0000 ; CALLING SEQUENCE: 
44 + 3 Square_root.wfc.v = MTHSCSORT (arg.rfc.r) 
9909 110 > INPUT PARAMETERS: 
0000 6 : The one input parameter is the address of a COMPLEX number (r, i) 
B60 BF : where r and i are both single-precision floating point values. 
0000 115 ; IMPLICIT INPUTS: 
SS 8 
0000 118 :; OUTPUT PARAMETERS: 
Bone Bs 3 NONE 
0000 121 : IMPLICIT OUTPUTS: 
oe igi 
900 124 ; COMPLETION CODES: 
44 : 5 : NONE 
9909 § > SIDE EFFECTS: 
000 ! : : Signals: Reserved Operand if r or i is bad (-0.0) 
00 131 :-- 
09 133 
000C 0 134 «ENTRY MTHSCSQRT, “M<R2, R3> 
bp 135 MTHSFLAG_JACKET ; flag as math routine 
6D 00000000'GF 9E 09 MOVAB G*MTHSSJACKET_HND, (FP) 
0 3 set handler address to jacket 
9 3; handler 
3 
52 Of BC 50 9 139 MOVE Bargadr (Ap). R2 3 RE zr 
52 8000 8F BA 0D 138 BICW #*x8000, R2 3; R2 = ABS(r) 


MTHSCSQRT 16-SEP-1984 01: AX/VMS Macro V04-00 Page 5 
1.008 MTHSCSQRT = compute COMPLEX square root pate} 9 Nf 38 YRTHRTL.§ SRCJ ATHCSORT.MAR: 1 . (4) 
mone Be ee a Wee FE niet em 
te _ MULF °5, RO : RO = (ABS(r) + CABS((r, 1))) / 2 
0990 "GF 1 1F 128 JSB G*MTHSSQRT_R3 ; RO = ROOT = SORT (aboves 
4 af p 3 14 MOVL argadr(AP), R2 3 R2 => (r, i) 
4 33 0058 148 event | ae. oe cheek 
1 04 D 146 CLRL eR : make zero quotient 
51 0602 50 4? O08) eb 18: BReRS AB. acre). Rt tare Ty ear 
oe ec en ne DE Ratt? * MO 
3 13) ‘ BGEQ RETRN > then ture’? (ROOT, Q) 
s Hone ie Ry Hae LTH lig inte 
Q 1 rH 13 BGEQ RETRNI > then Peturn <Q, ROOT) 
50 a. 4% 04 155 MNEGF R1, RO 3; else negate ROOT and Q 
51 53 $2 0047 136 MNEGF R3, R1 3; and return (-Q, -ROOT) 
0% foes 188 
0048 159 RETRN1: 
50 51 00 0048 #£«160 MOVL Ri, RO 3 continue to swap ROOT and Q 
51 53 DO OO4E 161 MOVL R3, R1 3 and return (Q, ROOT) 
ith 166 RETRN: 
04 51 16 RET 
005 164 
005 165 
005 166 END 
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ere = 00000004 

MTHSSJACKET_HND eeeeeee x 1 
MTHSCABS eerreere G 0 
MTHSCSQRT 00000000 RG 1 
mie 90000081 Root 
RETRNI 00000088 R 


+ 
! Psect synopsis ! 
Pees eeoeen ae se ase + 


PSECT name Allocation PSECT No. Attributes 
- (ABS. 00000000 ( 9. 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 00000052 (¢ 82.) G1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


Poemwooeecoeoceccccceccccn} 


Phase Page faults CPU Time Elapsed Time 
Initialization :00:00. :00:00. 
itializati 9 00:00:00.15 00:00:00.78 
pounens processing {4 00:00:00-49 Bp 80 B8-38 
Symbol table sort 2 80; Q: 8°83 00:00:09 :57 
Pass 2 4 00:00:00.49 8060: 0585 
Symbol table output Bh: 80 0075 0:00: “25 
Psect synopsis output 00:00:00.0 00:00:00.0 
Cross-reference output bs $4 88 88 09 «00 
Assembler run totals 300 00:00:02.09 0:00:11.49 


be working set Limit was 750 pages. 
2731 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages cf symbol table space allocated to hold 7 non-local and 2 local symbols. 
226 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
page of virtual memory was used to define 1 macro. 
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; Macro Library statistics ; 


Macro Library name Macros defined 
“S255$DUA28:(SYSLIBJSTARLET.MLB:2 0 (itt 0. 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE =SUPPRESSION/DJ SABLE=(GLOBAL , TRACEBACK) /LIS=LISS:MTHCSQRT/OBJ=OBJ$:MTHCSQRT MSRCS:MTHJACKE T/UPDATE=(ENHS :MTHJACKET) *MSRC 
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